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The Color-Flavor Transformation and Lattice QCD* 
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We present the color-flavor transformation for gauge group SU(A''c) and discuss its application to lattice QCD. 



1. Introduction 

Motivated by the study of disordered systems 
in condensed matter physics, the color-flavor 
transformation was first derived by Zirnbauer in 
1996 for an integral over U(7Vc) Q, 

In the above equation, -0 and -0 are Z2-graded 
tensors, i,j = l,...,Afc are color indices, a ~ 
1, . . . , n-i- and 6 = 1, . . . , ?^- are flavor indices, and 
the supermatrices Z and Z parameterize the coset 
space U(n+-|-n_|n+-|-n_)/U(n+|n+)xU(n_|n_). 
Their boson-boson parts satisfy Zbb = ^bb' 
and for the fermion-fermion parts, Zpp = — Zpp. 
The measure is given by D^a?^ {Z, Z) = d{Z, Z) 
Sdet(]L - ZZ)^". For further progress, see @. 

Notice that before the transformation, the color 
indices of tp and "ip are coupled by the matrix U, 
while the flavor indices are diagonal. In the trans- 
formed integral, the flavor indices are coupled by 
the matrix Z, and the color indices are diagonal 
— thus the name of the transformation. 

The integral over the gauge field closely re- 
sembles the one for a single link in the parti- 
tion function of lattice QCD at infinite coupling. 
We could therefore apply the transformation to 
all links of a lattice and then integrate out the 



fermions. In this regard, the most important 
property of the transformation ([l]) is that in terms 
of ip and -0, the right-hand side is local in co- 
ordinate space. Therefore, the fermion matrix, 
which for the usual lattice gauge action is an irre- 
ducible matrix of dimension ~ NcV, decomposes 
into blocks of size ~ Nf (number of flavors). This 
means that the fermion determinant of the trans- 
formed action is the product of determinants of 
fairly small matrices, which bears promise for a 
new approach to fermion algorithms. 

In order to apply this strategy to lattice QCD, 
we need to derive a color-flavor transformation 
for gauge group SU(A^c)- This is done in Sec. ||. 
In Sec. I we address the issue of how to go be- 
yond inflnite coupling (i.e. how to include the pla- 
quette action), discuss baryon loops that arise in 
SU(A'^c), and comment on simulation algorithms. 

2. Color-flavor transformation for SU(A^c) 

We consider two types of tensors '0 and ip with 
only fermionic degrees of freedom. This simplifles 
the calculations and is more appropriate for the 
application to lattice QCD that we have in mind. 
The reader may think of ip as corresponding to 
ipx+fL, "0 to <P to via;, and ip to ipx+ji- The 
result for the transformation for SU(Afc) is 



df/exp {i,lW^^i + ^lU^^^^^) 
dZdZ^ 



SU(JV, 



c^/x"/ det(l + ZZt)2A^/+^<= 
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*Based on two talks by the authors at Lattice 2002. 
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where a,b — 1, . . . , Nf, Xo — ^ 

XQ>o = Cq [det{Mf + det{Nf] . (3) 

The matrices in the determinants in Eq. carry 
color indices. They are M^^ = -0^(1 + ZZ'<)abipi 
and A/"'^ = (^^(1 + Z^Z)ab(pi. Finally, the con- 
stants are given by 

Co = Unio'^m + Nf + ny./{N,+ny.{Nj + ny. 
1 jt' {N, + ny.{Nf+ny. 



Q 



(4) 



The basic idea behind the proof of the color- 
flavor transformation is the following. One can 
define an auxiliary Fock space, with fermion cre- 
ation and annihilation operators acting in this 
space. The integrand of the SU(iVc) integral can 
be regarded as a complex-valued overlap of two 
states in this Fock space. Integration over SU(A^c) 
can be identified with a projection onto the color- 
neutral sector of the Fock space. The idea is to 
develop another implementation of this very pro- 
jector, which then leads to the color-flavor trans- 
formed integral. The details of the proof are 
rather complicated, and we shall only sketch some 
of the steps involved. For the full proof, we refer 
to Ref. |l (see also Ref. @). 

As mentioned above, we define two sets of 
fermionic creation and annihilation operators 
c^,c^ and da,^^ which carry color indices i = 
1, . . . ,Nc and flavor indices a = 1, . . . , Nf. These 
operators obey canonical anti-commutation rela- 
tions. The Fock space is obtained by acting on 
the Fock vacuum |0} with linear combinations of 
the creation operators. 

It is useful to define a certain "half-filled" state 
in this Fock space as 



IV'o) = (c}---4,)---(cf---c;^;)|o) . 

The outline of the proof is as follows. 



(5) 
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dC/exp {.piW^i^i + iplU^^^ipi) 
dU (V'o|exp(v;^< -c>'^ 



/SU(Afe) 



X exp 



{^o\exp{4>ldl,-cl^',) 

X Pexp (dlV'a+^aCa) IV^o) 



J2 ('Ao|exp(v;^< 



X 1q exp (^^4) IV'o) 



Cn 



D{Z, Zt) 



det(l -I- ZZt)^c 

Q=-Nf 

The interpretation of integration over SU(A^c) 
as a projection onto the color-neutral sector of 
Fock space can be understood by analogy with the 
case of projecting a vector in R'^ onto the z-axis. 
By averaging over all SO (2) rotations around the 
z-axis, what remains is the z-component of the 
initial vector, which is precisely the component 
that is invariant under such rotations. 

In order to derive an alternative expression for 
the projector P, we consider bilinear products 
of the creation and annihilation operators of the 
form Cj^c'g. Here, A,B — l,...,2Nf are com- 
posite indices introduced so that we don't have 
to distinguish between the c^'s and d^'s. These 
bilinears generate a g\{2NfNc) algebra, which 
has two commuting subalgebras that are of in- 
terest to us: the sl(iVc) color algebra and the 
gl(2iV/) fiavor algebra. The former is generated 
by operators = ^a'^a^A' ^ h '^'^'^ = 
Ea(caCa - c^I^c^^^), i = 1, . . . , iV, - 1, and the 
latter is generated by Eab = J^ii^ACg — ■^Sab)- 

States in the color-neutral sector satisfy 
T\Af) =0 for T G s\{Nc), or more exphcitly. 



Y:Z\S^a^aW) 



{Nf + Q) S'l \U) 



(7) 



The possible values for Q are Q = —Nf, . . . , Nf. 
Thus, the color-neutral sector splits into subsec- 
tors labelled by Q. Under the action of lJ{2Nf) 
on Fock space, given by 

g e Vi2Nf) ^ Tg^ exp(?^(logg)^BC*B) , (8) 

each sector Q subtends an irreducible representa- 
tion of the flavor group U(2iV/) with a rectangu- 
lar Young diagram that has Nf + Q rows and Nc 
columns. 
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Using the method of generahzed coherent 
states, we can express the projector onto the 
color-neutral sector as a sum of projectors onto 
the different subspaces Q, i.e. P = J2QL-Nf^Q■ 
To derive the projectors Ig, we start out with 
the highest weight vector of the representation 

Q, which is given by |V'q) = H^'^i UaL~[^ 
An overcomplete set of states is obtained by act- 
ing on this state with all group elements Tg, 5 G 
V{2Nf). In the sector Q = 0, the highest- 
weight vector I^/jq) corresponds to "half- filling." 
Clearly, Th |'0o) oc |^o) for h = diag(ft,+ , with 
h± € V{Nf ). The subgroup H = U(iV/) x U(iV/) 
of U(2Ai"/) is called the isotropy subgroup of |'0o)- 
In the sector Q = 0, coherent states may thus 
be parameterized without overcounting by the el- 
ements of the coset space \J{2Nf)/H. We label 
the elements of this coset space by picking a rep- 
resentative of each equivalence class gH, s(7r(g)). 
Each element g £ U(2A^/) can be decomposed 
into a product g — s{Tr{g))h{g). As an explicit 
expression we choose s(7r(g)) = 



{l + ZtZ)-l/2 -Zt(l + ZZt)-l/2 
Z(l ZtZ)-l/2 (l + ZZt)-l/2 



(9) 



where Z is an arbitrary Nf x Nf complex matrix. 
The action of Tg on \tpo) can be shown to give 

\Z) = cxpKZ,b4) \^o) det(l + 
where a,6 = 1, ... ,7V/. (10) 

The resolution of the identity in coherent states 
in the Q"Sector is given by 

1q = Cq dfiig) T, |^q) (V'qI . (11) 

In fact, this is a projector onto the Q-sector since 
it is easily seen to annihilate all other states. 

An integral over the whole group U(2iV/) 
can be decomposed into two integrals, one over 
H —\J(Nf)xlJ{Nf), and one over the coset space 
lJ{2Nf)/H, where the measure of integration 
over the coset space is given by D{Z,Z^) = 
CdZdZydet{l + ZZ^f^i. For Q = 0, the in- 
tegrand is indeed independent of h, and with the 
normalization yo\{H) = 1, the projector in this 
sector is simply 



The case of Q ^ is more complicated, since 
the integration over H is non-trivial. However, 
this integration can be done using standard group 
theory results, which finally leads to the color- 
flavor transformed integral of Eq. (||) . 

3. Application to Lattice QCD 

3.1. Induced QCD: How to generate the 
plaquette action 

The color-flavor transformation can be applied 
to each lattice link for lattice actions at infi- 
nite coupling, i.e. without a plaquette interaction 
term. Any program aimed at reformulating full 
QCD using the color-flavor transformation must 
therefore address the issue of including an inter- 
action term for the gauge fields. Here, we take ad- 
vantage of the idea of induced QCD, which goes 
back to Kazakov and Migdal Q and was further 
developed in Ref. |^ leading to the result of in- 
terest to us. The idea is to couple a number of 
additional heavy fermions to the gauge field. This 
induces gauge interactions, which can be directly 
related to an effective (non-zero) gauge coupling. 

To see how this works, we introduce Nh heavy 
fermions, described by the Wilson-Dirac operator 



±4 



5y^^+[,{r + "f^)Uf,{x) , (13) 



where we have not yet set r = 1. The hopping 
parameter k = \/{2Ma + 8r) tends to zero in 
the limit M ^ 00. We now integrate out the 
heavy quark fields to obtain the determinant of 
the Wilson-Dirac operator, D — 1 — kA. The re- 
sulting expression can then be expanded in pow- 
ers of the hopping parameter, 

dct^'-Z) = exp(iV,jtrlog(l - nA)) 



exp 



kA 



(14) 



Co / D{Z,Z^)\Z){Z\ 



(12) 



It is easily seen that terms containing odd pow- 
ers of A vanish when one takes the trace, since 
the matrix then has no entries which are di- 
agonal in lattice site indices. The quadratic 
term is just a constant (equal to zero for r ~ 
1), since U^ {x)Ul{x) = 1. Similarly, A'^ = 
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tr AxyAyzAzw^wx contains many constant terms 
as well as terms of the form 

tr [(r-7^)(r-7^)(r + 7^)(r + 7,J 

Ulixplix + v)U,{x + fL)U^{x)\ 



4(1 + - r*) tr Up 



(15) 



Here, p — (x; ^.v) denotes the lattice plaquette. 
Note the overall minus sign in the last line, which 
is essential to obtain the correct continuum limit. 
(In fact, with spinless fermions one obtains the 
wrong sign. The correct sign came from the lat- 
tice Lorentz structure of Wilson fermions, for any 
< r < 1. Staggered fermions also yield the cor- 
rect sign, which arises from the staggered phases. 
However, with staggered fermions, it turns out 
that one needs twice as many flavors to get the 
same coupling which correspondingly enlarges 
the dimension of the matrix Z in the color-flavor 
transformed action.) 

Collecting all relevant terms, we obtain at 
fourth order in k 



8Ar,,K*(l + 2r^ - r'^) V Re tr C/, 



(16) 



which is just the familiar plaquette action. We 
can thus identify 



(17) 



which becomes SNhK'^ for r = 1. 

In order to eliminate the effect of higher-order 
terms in the exponent at fixed lattice spacing a, 
one can let k ^ and — > oo in such a way 
that NhK^ = const. How large Nh has to be in 
practice to obtain the correct Yang-Mills theory 
must be determined numerically. A one loop cal- 
culation yields Nh > llNc/2 ||]. Incidentally, 
this implies that an approach based on the color- 
flavor transformation together with induced QCD 
is unlikely to be a suitable method for studying 
the large- iVc limit of Yang-Mills theory. 

3.2. Baryon loop expansion 

The color-flavor transformation deals with inte- 
grals that correspond to a single link of the lattice. 
In order to apply it to the lattice QCD partition 
function, we need to apply it to all links. Fur- 
thermore, the Dirac indices on -0 and ip have to be 



treated on par with the flavor indices as far as the 
transformation is concerned. Thus, the dimen- 
sion of the Z matrices obtained after the transfor- 
mation is given by dim(Z) = 4(7V/ -|- Nh) = 47V,. 

For Wilson fermions, the action at inflnite cou- 
pling is given by 



s 



2 ^ 

x,fi>0 



x.a 



(m,-f4r)4";>"': 



(18) 



We have included all the indices, namely Greek 
Dirac indices a, P = 1, ... ,4, color indices i,j = 
1, . . . ,Nc and flavor indices a = 1, . . . , Nq. Note 
that the 7^ are traceless, idempotent matrices, 
which can be diagonalized by unitary matrices, 

7;. = r^diag(i,i,-i,-i)rt . 

If we now set r = 1, we have 
1 



(19) 



{r + 7m) 
{r - 7m) 



r p, pt 



(20) 

(21) 



Here, P+ = diag(l, 1, 0, 0) and P_ =diag(0, 0, 1, 1) 
are projection operators in Dirac space. On each 
link, we define 



Yx,a^ /J. ^ — 



'Px,a — + ^ ti 'Px,a ; Vx-^ 



Vx+il,a^ M + 



t;+M,a — M ^x-\ 

This definition is of course not the only possible 
choice for the inclusion of the Dirac matrices, and 
therefore the color-flavor transformed action will 
not be unique. This should, however, have no 
influence on the physics, since the transformation 
is in any case still exact. The (f) and (f) above 
are still Grassmann variables, and so we can now 
apply the color-flavor transformation to each link 
of the lattice to obtain a transformed partition 
function of the form 

dZ^{x)dZl{x) 



^QCD 



n 



dct{l + Z^{x)zl{x)) 



-s^ix) 

Q=0 



XQmW- (22) 
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In the exponential, there is now a local term 
S^{x) = ^^:Mxrjr,';, where 



(23) 



Furthermore, we have XOfii^) = 1 ^-nd 
XQ^(x)=CQ[det'5A^^(x)+det«A/;(a:)] . (24) 
The matrices which appear here are defined as 
M^{x) = 

^x,a^ p ^- \^ ^ '^x.4i^x,p) ^- ^ p 'yx+fi,b ■ 

The next step is to integrate out the fermions. 
This is straightforward for Q — 0. With the def- 
inition of B^{x) from above, we define B{x) ~ 
J2 Bp{x) and B = B{x). Then B has a block- 

/I all sites X 

diagonal structure, with different blocks corre- 
sponding to different lattice sites. After integra- 
tion over the Grassmann variables, we thus obtain 



7Q=o 

^QCD 



det^^B = Yl det^'^B(x) 



(25) 



all sites X 



For Q — 0, the situation becomes more com- 
plicated, since the XQ>Q terms induce non-local 
contributions to the action. Notice that these 
terms correspond to baryon propagation, which 
explains why they were absent in the case of gauge 
group \J{Nc). 

The case of one static baryon (Q = 1) was stud- 
ied in saddle-point approximation (for large Nc) 
in 1-1-1 dimensions by Budczics et al. in Ref. [Q. 
However, our goal is an exact (numerical) ap- 
proach to QCD, which takes into account all pos- 
sible values of Q. 

In order to evaluate terms from XQ>Q contri- 
butions, we need to apply Wick's theorem. This 
results in terms involving the inverse of the ma- 
trix B. Note that B is diagonal in space and 
color, so that we have 



(26) 



= 1 



Figure 1. Expansion of baryon contributions: 
Only closed loops survive the integration over the 
Grassmann variables. 

where p^q — 1, . . . , 4:Nq combine flavor and Dirac 
indices. Let us denote 



x + p. X x + p. 



(27) 



Here, a solid link denotes the presence of 
1,2,..., ANq baryons. Now consider 



n 



1 



X + fl 



X + fl 



(28) 



We can expand the product and note that be- 
cause of {B~^)xy ^ Sxy, only those terms survive 
that correspond to closed loops (cf. for example 
Fig. Hence, in the partition function only 
closed baryon loops are allowed. These must not 
backtrack, so that each link occurs at most once 
in a graph. The paths may self intersect, and dis- 
connected loops are also allowed. In addition, for 
each site the total number of baryons on the links 
connected to this site must not exceed 4A^g. 

Let us now consider observables, such as cor- 
relation functions. (In the color-flavor trans- 
formed formulation, gauge fields no longer ap- 
pear, and it is therefore not surprising that non- 
gauge-invariant quantities are manifestly equal to 
zero.) As an example, take the baryon two-point 
function (ByBz), with B = Sijkip^ip-'ip^ ■ Schemat- 
ically, we have to compute 



a:, /J, 



1 



-t- 



X + jl X 



X + fl 



(29) 



Contributing baryon paths are those which orig- 
inate at the site y and terminate at site z. Of 
course, a background of closed baryon loops, such 
as those appearing for the partition function, is 
also allowed, in addition to the baryon propagat- 
ing from y to z. 
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3.3. Algorithmic considerations 

While it is, in principle, clear how to perform 
a simulation for the Q = case, the presence of 
baryon loops poses additional problems. Clearly, 
the number of possible loops grows prohibitively 
with the lattice volume so that a probabilistic ap- 
proach to the generation of such loops is needed. 
A possibility is some form of random walk al- 
gorithm 0. It should be noted, though, that 
the coefficients Cq decrease rapidly as functions 
of Q so that multi-baryon paths are potentially 
strongly suppressed. For example, for Q = 1 and 
A^c = 3, we have Ci cx l/N^, which is already 
much smaller than in the Q = case. In addition, 
since longer loops involve more factors of Cq, only 
small loops, or shortest paths in the case of corre- 
lation functions, will significantly contribute. An 
expansion in baryonic loops may therefore be a far 
more systematic and better controlled approach 
than quenching, for example. 

Unfortunately, the color-flavor transformed ac- 
tion is complex. This is already manifest for Q = 
0, as can be seen from Eqs. (p^ ) and ([2^), since 
Z is a general complex matrix. We have tried to 
attack this problem by partially integrating out 
some of the degrees of freedom of Z (cf. Ref. [D), 
and also by finding different parameterizations of 
the coset space U(2iV/)/[U(7V/)xU(iV/)], but so 
far have not been able to obtain a satisfactory 
result. However, some progress has been made 
recently in attacking fermion sign and complex 
action problems |^, and this may also turn out 
to be helpful for the color-ffavor transformation. 

4. Summary and Outlook 

We have derived the color-flavor transforma- 
tion for gauge group SU(A^c)- This transfor- 
mation converts a certain type of integral over 
color gauge flelds to an integral over flavor ma- 
trices. The integral being transformed is of the 
same form as the one encountered on each link 
of the lattice QCD partition function at infinite 
coupling. Since (for zero baryons) the fermion 
determinant after the transformation is block- 
diagonal in the lattice site indices, this approach 
has the promise of leading to new ways of sim- 
ulating dynamical fermions. However, in con- 



trast to gauge group U(7Vc), the SU(A^c) case 
gives rise to additional baryon loops, which lead 
to nearest-neighbor coupled fermion terms. Nev- 
ertheless, multi-baryon loops are potentially sup- 
pressed, and a systematic expansion of the parti- 
tion function in terms of such loops is possible. 

In order to go beyond the infinite coupling 
limit, a number Nh of heavy flavors can be added 
to the theory to generate the plaquette action 
{g^^ — SiV/jK^ for Wilson fermions). 

In the present formulation, the color-flavor 
transformed action is complex, which makes nu- 
merical simulations based on this approach im- 
practical without further improvements. Work in 
this direction is in progress. 
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